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Abstract 

We present a measurement of the branching fraction for the flavor-changing neutral current 
process B — > Xg £^£~ based on a sample of 88.9 x 10^ T(4S') BB events recorded with the 
BABAR detector at the PEP-II e~^e~ storage ring. The final state is reconstructed from pairs of 
electrons or muons and a hadronic system consisting of one or Kg and up to two pions, with 
at most one vr*^. Summing over both lepton flavors, we observe a signal of 41 it W{stat) =t 2{syst) 
events with a statistical significance of 4.6 a. The inclusive branching fraction is determined to be 
B{B Xs e^e-) = (6.3 ± l.6{stat)+\i{syst)'^ x lO'^ for m(£+^-) > 0.2 GeV/c^. Ah results are 
preliminary. 
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1 Introduction 



The rare decay B — s- Xg proceeds through a b ^ si~^i~ transition, which is forbidden at 

tree level in the Standard Model (SM). However, such a flavor-changing neutral current (FCNC) 
process can occur at higher order via electroweak penguin and W~^W~' box diagrams, as shown in 
Figure n The h s£'^i~ transition therefore allows deeper insight into the effective Hamiltonian 
describing FCNC processes and is sensitive to the effects of non-SM physics that may enter these 
loops; see, for example, Refs. Hill]. 

Recent calculations of the branching fractions for the exclusive decay modes predict B{B 
Ke+e-) = B{B Kfi+fi-) = (0.35 ± 0.12) x lO^^, B{B K*e+e-) = (1.58 ± 0.49) x lO"*^ 
and B{B K*fi+fi-) = (1.19 ± 0.39) x 10"^, and for the inclusive processes B{B Xs e+e') = 
(6.9±1.0) X 10"*^ and B{B X^/i+zx") = (4.2±0.7) x 10"^ p. In the electron channel, the branch- 
ing fraction is predicted to be B{B Xs e+e") = (4.2 ±0.7) x 10"^ |2l for m(e+e-) > 0.2 GeV/c^. 
Measurements of the inclusive branching fractions are motivated by the smaller theoretical uncer- 
tainties, as compared with predictions for exclusive decays. Both the Belle and BABAR Collabo- 
rations have observed the B — > Ki~^i~ {I = e,^) decay jUIHl, BABAR has found evidence for the 
B — > K*£~^i~ decay |6i, and Belle has measured the inclusive B ^ Xs £~^£~ decay Q. 

In the present analysis, we study the B Xsi~^i~ process by reconstructing the final state from 
pairs of electrons or muons and a hadronic system consisting of one or and up to two pions, 
with at most one vr'^. The choice of maximum number of pions is motivated by the fact that the 
signal-to-background ratio drops significantly with increasing multiplicity. This approach [S] , which 
sums over exclusive modes, allows approximately half of the full inclusive rate to be reconstructed. 
If the fraction of modes containing a is taken to be equal to that containing a Kg, the missing 
states that remain unaccounted for represent ~25% of the total rate. We require the hadronic mass 
for the selected final states to be less than 1.8 GeV/c^ to reduce combinatorial background. This 
cut retains approximately 94% of the signal in the reconstructed modes. We correct for the missing 
modes and the effect of the hadronic mass cut to extract the inclusive B ^ Xs I'^i" decay rate for 
m{£+£^) > 0.2 GeV/c2. 

b 1 ' » s b 

u,c,t 

Figure 1: Feynman diagrams for h — > sl'^l~ transitions. 



2 The BaBAR Detector and Data Sample 

The data sample used in this analysis was collected with the BABAR detector during 1999-2002 at 
the PEP-II asymmetric-energy e^e~ storage ring at the Stanford Linear Accelerator Center. It 
consists of 81.9 fb^^ recorded at the T(45') resonance, corresponding to 88.9 x 10^ 7(45) BB 
events, as well as an additional 9.6 fb"^ recorded 40 MeV below the r(45) resonance. The off- 
resonance data provide a clean sample of e~ — > qq (with q = u,d,s,c) continuum events for 



W"/ 

L 



r r 



u,c,t 
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background studies. 

This analysis relies on the charged-particle tracking systems and particle identification systems 
of the BABAR detector j^j. A 5-layer silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) provide tracking and particle identification for charged particles. The DIRC, a Cherenkov 
ring-imaging system, is primarily used for charged hadron {7r,K,p) identification. The electro- 
magnetic calorimeter (EMC), consisting of Tl-doped Csl crystals, provides electron identification 
and photon reconstruction. These detectors are located inside a 1.5-T solenoidal superconducting 
magnet. Muons are identified in the instrumented flux return (IFR) with resistive plate chambers 
interleaved with iron plates. 

Clean particle identiflcation for the final state particles e^, /x^, K^, Kg, tt^ and ttO is important 
for this analysis. Electrons and muons are required to have lab-frame momenta greater than 
0.5 GeV/c and 1.0 GeV/c, respectively. Bremsstrahlung photons from electrons are recovered 
by combining an electron with up to three photons within a small angular region around the 
electron direction ^U]. candidates are reconstructed from pairs of oppositely-charged tracks 
with |m(7r^7r^) — m{Kg)\ < 11.2 MeV/c^, decay length greater than 2 mm, and angle between 
the Kg momentum vector and the line between the primary vertex and the Kg vertex satisfying 
cos (5 > 0.99. Charged pions are selected from tracks after rejection of cleanly identified and K^ 
candidates. Neutral pions are required to have lab- frame energy greater than 400 MeV, photon 
daughter energies greater than 50 MeV, and a 77 invariant mass that satisfies |m(77) — m{'ir^)\ < 
10 MeV/c2. 

3 Analysis overview 

We use a technique that sums over exclusive modes to measure the inclusive branching fraction 
B{B Xg Compared to a fully inclusive approach, this method has the advantage of 

being able to exploit the strong kinematical discrimination provided by m^s = E'^eam ~ Pb 
AE = Eb — Ebeam, where Ebeam is the beam energy and Eb (pb) is the reconstructed B meson 
energy (3-momentum) . All quantities are evaluated in the e^e~ center-of-mass system (CM). This 
analysis approach is needed to suppress large backgrounds from BB and continuum events, and to 
extract the expected small signal yield. 

The main contribution to the combinatorial background is from semileptonic decays in BB 
events. In these events, B Xgi~^i~ candidates contain decay products from both B mesons (i.e., 
the daughters of the B ^ Xg i'^l^ candidate do not originate from the same point) and there is a 
significant amount of missing energy due to neutrinos. Another contribution to the combinatorial 
background is due to continuum events, which are effectively suppressed with event-shape variables. 

The dominant peaking backgrounds, which are from B — > J/tpX and B — > ip{2S)X decays with 
J/'ip{ip{2S)) i~^£~ , mimic the signal and need to be efficiently removed with cuts on the dilepton 
mass m{£~^i~). The resulting veto sample provides a large control sample of decays with a signature 
identical to that of the signal, albeit in a restricted range of m{e+e-). B L»(*)n7r (n > 0) decays 
with misidentification of two charged pions as leptons present an additional source of peaking 
background. 

Simulated B ^ Xg £'^i~ events are produced with a combination of exclusive and inclusive 
models. In the hadronic mass region of m{Xg) < 1.1 GeV/c2, exclusive B K^*k+i- decays 
are generated according to ^ , where the relevant form factors are computed using light-cone 
QCD sum rules. The remaining decays, in the region 'm{Xg) > 1.1 GeV/c^, are generated with 
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a nonresonant model following ^ and using the Fermi motion model of 0^ • JETSET is 
then used to hadronize the system consisting of a strange quark and a spectator quark. 

4 Event selection 

Events are required to have a well determined primary vertex, be tagged as multi-hadron events, 
have a ratio of the second to zeroth-order Fox- Wolfram moments ^Hl (calculated using charged 
tracks and neutral calorimeter clusters) of R2 < 0.5, have at least four charged tracks, and contain 
at least two leptons (either electrons or muons) satisfying loose particle identification criteria and 
having lab-frame momenta greater than 0.5 GeV/c for electrons and 0.8 GeV/c for muons. 

In the next step, B candidates are reconstructed. A dilepton is selected by choosing the e+e~ 
or pair with the largest value of |p(^'^)| + |p(£^)|, using lab-frame momenta. The dilepton 

candidate is required to form a good vertex. Using this (.^i" pair, B Xsi~^i~ candidates are 
formed by adding either a or a and up to two pions, but no more than one tt^. In this 
manner, ten different hadronic topologies are considered: K^, K~^t:^, K~^'k~ , K~^Tr~Tr^, K~^Tr~Tr~^, 
K% K^TT°, K^TT+, K^TT+TT°, and K^tt+tt-. 

The resulting B candidates are required to have an vertex fit probability log{Puytx) > —10, 
a B vertex constructed from charged daughter particles with fit probability log{PBvtx) > —10, 
m{Xs) < 2.5 GeV/c^, 5.0 < mss < 5.29 GeV/c^, and \AE\ < 0.3 GeV. 

At this stage, there is an average of five B candidates per event in the signal simulation. Only 
the B candidate with the largest signal likelihood value is retained. The signal likelihood function 
is constructed based on the simulated distributions of AE, log{PBvtx), and cosOb for the signal, 
where is the angle between the momentum vector of the B candidate and the beam axis in the 
CM frame. We select the best candidate before further selection cuts are applied. This approach 
entails a small cost in signal efficiency, which is more than compensated by an improved background 
suppression in the final sample. 

After selecting the best candidate, combinatorial backgrounds are suppressed by tightening 
some of the earlier requirements: m{Xs) < 1.8 GeV/c^, 5.20 < tues < 5.29 GeV/c^, and —0.2 < 
AE < 0.1 GeV. In addition, we require the separation between the two leptons along the beam 
direction to satisfy \Az\ < 0.15 cm, where the z-coordinate of each lepton is determined at the 
point of closest approach to the beam axis. We also require the rest of the event, formed by 
combining all charged tracks and neutral calorimeter clusters not included in the B candidate, to 
satisfy -5.0 < AE^'^^ < 2.0 GeV and > 4.9 GeV/c^. 

Charmonium backgrounds are reduced by removing B candidates with dilepton mass in the 
ranges 2.70 < m(e+e-) < 3.25 GeV/c^, 2.80 < m(/i+/x-) < 3.20 GeV/c^, 3.45 < m{e+e-) < 
3.80 GeV/c^, and 3.55 < m{n+ < 3.80 GeV/c^. If one of the electrons from a J/ip or V(25) 
decay erroneously picks up a random photon in the Bremsstrahlung-recovery process, the dilepton 
mass can increase sufficiently to evade the above cuts. Therefore the charmonium veto is applied 
to the dilepton mass before and after Bremsstrahlung recovery. Using the simulation, we estimate 
the remaining peaking charmonium background to be 0.6 it 0.3 events and 0.3 ±0.2 events for e~^e~ 
modes and modes, respectively. 

The potential peaking background from B — > X^^ decays, followed by conversion of the photon 
into an e^e~ pair in the detector material, is a concern for the e^e^ modes only. After observing an 
excess of events in data over the simulation for low e^e~ mass in the tues < 5.27 GeV/c^ sideband 
distribution, we completely remove this background by requiring m{e~^e'') > 0.2 GeV/c^. 

The final suppression of the combinatorial background is achieved with a likelihood ratio 
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based on nine variables: AE, AE^'^^ , m^g^, Az, log{PBvtx), 

cos O^igg , where O^^gg is tlie angle 
between the missing momentum vector for the whole event and the z axis in the CM frame, cos Ob, 
I cos I) where 9t is the angle between the thrust axes of the B candidate and the rest of the 
event in the CM frame, and R2. The variables AE, AE^^^ , and fn^^ are most effective at 
rejecting BB background, especially for events with two semileptonic decays, which have large 
missing energy. For continuum suppression, the event-shape variables |cos^r| and R2 are most 
useful. The distribution for the background-discrimination variable i is parameterized by a PDF 
rH for each component j of the sample: j = 1 for signal, j = 2 for BB background, and j = 3 for 
continuum background. A likelihood is then computed for each sample component as the product 
of the individual PDFs for that component: = IliVi- The ratio between the likelihood for the 
signal component and the sum of the signal and background likelihoods then provides the final 
discriminating variable, the likehhood ratio £r = c^ignai / f^^signai _^ ^bb _^ £cont^^_ 

Using the simulation, the cut on Cji is optimized to maximize the statistical significance of the 
signal. This optimization is performed separately for e~^e~ modes and modes in the three 

m{Xg) regions m{Xg) < 0.6 GeV/c^, 0.6 < m{Xg) < 1.1 GeV/c^, and 1.1 < m{Xg) < 1.8 GeV/c^, 
resulting in the cuts Cr > 0.3, 0.4, and 0.9 for the e+e~ modes and Cr > 0.2, 0.6, and 0.9 for the 
fi'^fi~ modes. 

After applying all selection criteria, a sample of 349 B ^ Xg e^e^ and 222 B — > Xg fJ,^fJ-~ 
candidates remains. According to the simulation, the background left at this stage of the analysis 
consists mostly of BB events (86% and 72% of the total background in the electron and muon 
channels, respectively). 



5 Maximum likelihood fit 

We perform an extended, unbinned maximum likelihood fit to the tues distribution in the region 
iT^ES > 5.2 GeV/c^ to extract the signal yield as well as the shape and yield of the combinatorial 
background. The likelihood function C consists of four components: 

C = n ii^srg + N,-peak)n-' + N.-peakV^'''''' + N^ombV^'] 

k=l 

(1) 

where and k denote the total number and index of candidate events, respectively. Ngig, N^-peak, 
Xh-peak-, and Ncomb represent the yields of the signal, charmonium peaking background, hadronic 
peaking background, and combinatorial background, respectively, with corresponding PDFs given 
by -p'^P^°'^ ^ -ph-peak ^ jycomb ^ The signal PDF is the same as that for the charmonium peaking 
background since it is extracted from the charmonium-veto data sample. 

The signal PDF "P^ ^'^"'^ is described by a Gaussian for /^^/i^ modes as well as for e^e~ modes, 
since the Bremsstrahlung recovery and selection procedure for e^e~ modes lead to a negligible 
radiative tail in the tues distribution. The Gaussian shape parameters are determined from fits 
to the sum of a Gaussian and Argus JHI function for the charmonium-veto data sample. The 
fits yield signal mEs peak positions and resolutions of nigig = 5280.04 ± 0.05 MeV/c^ and agig = 
2.80 ± 0.05 MeV/c^ for the e+e" modes, and rugig = 5280.05 ± 0.07 MeV/c^ and agig = 2.61 ± 
0.06 MeV/c^ for the fJ'^fJ'~ modes, respectively. In the simulation, the Gaussian fit results for the 
niES distributions for correctly reconstructed signal are in agreement with the shape parameters 
extracted from the fits to the charmonium-veto sample. The signal yield Ngig is a free parameter 
in the likelihood fit. 
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The charmonium peaking background is estimated from the simulation to contribute N(.-peak = 

0.6 lb 0.3 events in the electron modes and N(.-peak = 0.3 it 0.2 events in the muon modes. In the 

likelihood fit, Nc~peak is fixed to these values. 

The size and shape of the hadronic peaking BB background component arising from B — >■ 
rnr (n > 0) decays with misidentification of two charged pions as leptons are derived directly 

from data by performing the analysis without the lepton identification requirements. Parameters 

for the PDF p^-P^'^'^ 

are taken from a fit to the sum of a Gaussian and Argus function, with uies 
peak positions and resolutions of mp^ak = 5280.13 ± 0.06 MeV/c^ and cjpf.ak = 2.99 ± 0.07 MeV/c^. 
Taking the vr to ^ misidentification rates into account, the remaining hadronic peaking backgrounds 
are estimated to be A'^fj.peafe = 2.4 it 0.8 events for the fJ'~^fJ^~ modes and are negligible for the e~^e~ 
modes. In the likelihood fit, Ny^^p^ak is fixed to the estimated values. 

The combinatorial background PDF, 'p^"^^^ is given by an Argus shape, which describes the 
combinatorial contribution from continuum events and BB events. The Argus cutoff is determined 
by the beam energy in the T(4S') rest frame, Ebeam = 5.290 GeV, whereas the Argus shape 
parameter and the yield Ncomb are free parameters in the likelihood fit. 



6 Results 

Using the fit parameterizations described above, we fit the rriES distributions for the selected B 
Xse~^e~ and B candidates separately and obtain the results shown in Figured The fit 

results are summarized in Tabled The statistical significance is 5 = 2{ln Cmax — ^^J~-max)j where 
^max represents the maximum likelihood for the fit and C^^^ denotes the maximum likelihood for a 
different fit when the signal yield is fixed at Ngig = 0. The significance does not incorporate effects 
due to systematic uncertainties in the signal shape or in the contribution from peaking background. 
However, these have been found to be small. The B — > Xg £~^£~ signal yield presented in Tabled 
is the sum of the B — > Xs e'^e~ and B ^ Xg l-i'^ fJ-^ signal yields. A separate fit to the combined 
electron and muon channels gives a very similar result. 

Figure |2{d) shows the distribution of iriES for B ^ Xs e^/x^ candidates selected using the 
nominal selection criteria but requiring that the two leptons have different flavor. This sample 
provides a cross-check for the parameterization of the combinatorial background. Figure Efd) 
shows an Argus fit to the niEs distribution and, as expected, there is no evidence for peaking 
background. The Argus shape and yield parameters for the combinatorial background are free in 
the nominal fit to data. 

Figure |31 shows the distribution of m{Xs) for electron and muon channels combined. This 
distribution is obtained by performing the nominal likelihood fit in separate m{Xs) regions. Figure|31 
indicates that the observed signal receives contributions from final states across a range of hadronic 
mass, including hadronic systems with mass above that of the i^*(892). 

The branching fraction B for the signal is calculated from 

where iV^^ = (88.9 ± 1.0) x 10^ is the number of BB pairs produced in 81.9 fb ^ and e is the signal 
efficiency. The fitted signal yield Ngig contains a contribution from misreconstructed B — > Xg £^i~ 
decays (cross-feed events) of l.Sitl.S events in the electron channel and 0.9ib0.9 events in the muon 
channel. The central values are derived from fits to data in which the nominal fit function is changed 
to include a cross-feed PDF whose shape and normalization with respect to that for the correctly 
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m^s ( GeV/c^ ) m^g ( GeV/c^ ) 



Figure 2: Distributions of mEs for selected (a) B ^ Xg e+e~, (b) B — > Xg fi^fj.^ , (c) B — > Xsl^i~ 
{(. = e,id), and (d) B — > e^fi^ candidates. The solid lines represent the result of the fits and 
the dashed lines represent the background component under the signal peaks. 



reconstructed B Xgi'^i^ decays are taken from the simulation. The corresponding uncertainties 
are estimated from Monte Carlo studies using a large number of data-sized experiments. The signal 
efficiency e is adjusted to reflect the contribution from cross-feed events. 

7 Systematic uncertainties 

Systematic uncertainties are of two different types: those that affect the extraction of the number 
of signal decays and those that affect the calculation of the branching fraction. The systematic 
uncertainties are summarized in Table |51 

Uncertainties affecting the extraction of the signal yield are evaluated by varying the signal 
Gaussian parameters (mean and width) and the signal shape (using asymmetric signal shapes) 
within the constraints allowed by the charmonium veto data. The small amount of peaking back- 
ground from B — >■ J /ip X and B ip{'^S) X decays is varied according to the estimated number of 
such decays: 0.6 it 0.3 and 0.3 ± 0.2 in the electron and muon channels, respectively. The peaking 
background in the muon channel due to B ^ Xs vr^vr" decays is varied according to the estimated 
number of such events: 2.4 it 0.8, where the uncertainty is dominated by the uncertainty in the 
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Table 1: Results of the fit to data: signal yield, charmonium and hadronic peaking backgrounds 
(fixed in the fit), combinatorial background yield, and statistical significance. 



Sample 


Ns^g 


Xc—peak 


X h—peak 


Xcomb 


Signif. 




29.0 ± 8.3 


0.6 


0.0 


319.4 ± 18.9 


4.0 


Xs iJi'^l-r 


12.4 ± 6.2 


0.3 


2.4 


207.0 ± 15.2 


2.2 


Xs £+r 


41.4 ± 10.3 


0.9 


2.4 


526.4 zb 24.3 


4.6 



!5 



25 



f 20 



15 
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Figure 3: Distribution of number of signal events as a function of hadronic mass for electron and 
muon channels combined for data (points) and Monte Carlo signal (histogram). The vertical error 
bars represent statistical uncertainties only. 



vr to /i misidentification rates determined from control samples in the data. The shape of this 
background is also varied according to the data sample selected without lepton identification. The 
largest variations in the signal yield are observed for changes in the signal shape and in the amount 
of peaking background. 

Uncertainties affecting the signal efficiency originate from the detector modeling, from the 
simulation of signal decays, and from the estimate of the number of B mesons in the sample. By 
far the dominant component is that due to the simulation of signal decays, discussed in detail 
below. 

The detector modeling uncertainty is sensitive to the following uncertainties determined from 
the data: the uncertainty in the tracking efficiency of 1.3% per track (2.0% per charged pion); the 
uncertainty in the charged-particle identification efficiency of 0.7% per electron, 1.6% per muon, 
1.0% per kaon and 2.0% per pion; and the uncertainty in the reconstruction efficiency of 10% per 
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Kg and 5% per ir^. The efficiency of the Ukehhood ratio cut to suppress combinatorial background 
is checked with the charmonium-veto sample and the level of discrepancy with the simulation is 
taken as the corresponding uncertainty. The fraction of signal cross-feed included in the signal 
Gaussian is studied with a parameterized Monte Carlo simulation technique and the uncertainty 
is set to correspond to the R.M.S. obtained for a large number of data-sized experiments (the 
uncertainty amounts to a variation of ±100% of the number of signal cross-feed events included in 
the Gaussian). 

The dominant source of uncertainty arises from modeling signal decays. Parameters of the Fermi 
motion model are varied in accordance with measurements of hadronic moments in semileptonic B 
decays |17l and the photon spectrum in inclusive B — > Xg 7 decays ^HI- The fraction of exclusive 
B Ki'^l^ and B K*l~^l~ decays is varied according to theoretical uncertainties and the 
corresponding systematic uncertainties are added linearly (i.e., we assume the uncertainties in the 
branching fraction for these two modes are fully correlated). The transition point in m{Xs) between 
pure K*i~^l~ and nonresonant Xg £~^i~ final states is varied by ±0.1 GeV/c^. 

The nonresonant Monte Carlo event generator relies on JETSET to fragment and hadronize 
the system consisting of a final state s quark and a spectator quark from the B meson. Since the 
signal efficiencies depend strongly on the particle content of the final state, uncertainties in the 
number of charged and neutral pions and in the number of charged and neutral kaons translate 
into a significant uncertainty in the signal efficiency (for m{Xs) > 1.1 GeV/c^). 

The ratio between the generator yield for decay modes containing a Kg and that for modes 
containing a charged kaon is varied according to 0.50 ± 0.05, to allow for isospin violation in the 
decay chain. The ratio between the generator yield for decay modes containing one vr*' and that 
for modes containing none is varied according to 1.0 ± 0.5. The ratio between the generator yield 
for decay modes into two-body hadronic systems and that into three-body hadronic systems is 
varied according to 0.5 ±0.3. Uncertainties in the last two ratios are set by the level of discrepancy 
between data and Monte Carlo as measured in the BABAR semi-inclusive B — > Xg 7 analysis |19j . 
(Both B Xg 7 and present analyses use the same hadronization model.) 

The ten modes selected in this analysis only capture about 50% of the full set of final states. 
Approximately half of the missing modes is due to final states with a K^ meson and their contri- 
bution can be determined from the Kg modes. However, we need to account for the uncertainty 
in the fraction of modes with too many pions or kaons (two extra kaons may be produced via ss 
popping), as well as for modes with photons that do not originate from vr*^ decays but rather from 
r/, rj' , Lj, etc. For final states with m{Xs) > 1.1 GeV/c^, we vary these different fractions by ±50%. 

Including systematic uncertainties, the measured branching fractions for m{£~^i~) > 0.2 GeV/c^ 

are 



where the first error is statistical and the second error is systematic. The combined B Xg i~^£~ 
branching fraction is the weighted average of the branching fractions for the electron and muon chan- 
nels, where we assume the individual branching fractions to be equal for m{£^i^) > 0.2 GeV/c^. 
Table |21 summarizes the results of the analysis and lists both the statistical and systematic errors 
in the signal yields, the signal efficiencies and the branching fractions. 




(3) 
(4) 
(5) 
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Table 2: Summary of fractional systematic uncertainties (in %). The uncertainties in extracting the 
signal are presented first and those related to the signal efficiency and BB counting are presented 
second. 



bource 


Xs e^e 


V + — 
Xs iTjl 


Signal shape 


±4.3 


±2.5 


Peaking background 


±1.0 


±5.8 


Signal yield total 


±4.4 


±6.3 


Tracking efficiency 


±5.0 


±4.8 


Lepton identification efficiency 


±1.3 


±3.2 


Kaon identification efficiency 


±0.8 


±0.8 


vr^ identification efficiency 


±1.5 


±1.3 


Kg efficiency 


±2.2 


±1.9 


7r° efficiency 


±0.8 


±0.5 


LR cut efficiency 


±3.7 


±3.3 


Cross-feed efficiency 


±6.7 


±7.4 


Detector model subtotal 


±9.7 


±10.2 


Fermi motion model 


4-iU.2 
-4.0 


+iU.8 
-3.6 


B{B K£+e~) 


±9.6 


±12.6 


B{B K*£+e-) 


±6.2 


±6.8 


K*-Xs transition 


±3.3 


±4.1 


Hadronization 


±11.1 


±7.5 


Missing modes 


±10.6 


±8.9 


Signal model subtotal 


+24.5 
-22.6 


+25.4 
-23.3 


Monte Carlo statistics 


±2.5 


±3.0 


BB counting 


±1.1 


±1.1 


Efficiency and total 


+26.5 
-24.8 


+2V.6 
-25.6 



8 Summary 

Using a sample of 88.9 x 10^ r(4S') BB events, we measure the branching fraction for the rare 
decay B — > Xs i^i^ , where £ = e or /x, with a sum over ten different hadronic states (with up to 
two pions). For m{£~^£~) > 0.2 GeV/c^, we observe a signal of 41 ± W{stat) ± 2(syst) events and 
obtain a branching fraction of 

B{B Xs i+r) = (6.3 ± 1.6{stat)t\i{syst)'^ x 10"^ (6) 

with a statistical significance of 4.6 a. This preliminary measurement agrees well with the recent 
prediction by Ali et al. ^ and the measurement performed by the Belle Collaboration, B{B — > 
Xs = (6.1 ± lA{stat) ±li (syst)) x lO"*' for m{£+£-) > 0.2 GeV/c^ p. 
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Table 3: Summary of results: signal yield, statistical significance, efficiency and branching fraction. 
In the case of the signal yield and the branching fraction, the first error is statistical and the second 
error is systematic. In the case of the signal efficiency, the first error corresponds to uncertainties in 
detector modeling, BB counting, and Monte Carlo statistics, whereas the second error corresponds 
to the uncertainties in the signal model. 



Sample 


N ■ 
sig 


Signif. 




B (xlO-*^) 


Xs e+e- 
Xg jj+jj.- 
Xs 


29.0 ± 8.3 ±1.3 
12.4 ± 6.2 ±0.8 
41.4 ± 10.3 ± 1.5 


4.0 
2.2 
4.6 


2.46 ± 0.25^^-^^ 
1.23±0.13t°J^ 


6.6±1.9;|:^ 
5.7±2.8;j;^ 
6.3±1.6li;^ 
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